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We discuss the primordial spectrum of a massless and minimally coupled scalar field, produced 
during the initial anisotropic epoch before the onset of inflation. We consider two models of the 
anisotropic cosmology, the (planar) Kasner de Sitter solution (Bianchi I) and the Taub-NUT de 
Sitter solution (Bianchi IX), where the 3-space geometry is initially anisotropic, followed by the de 
Sitter phase due to the presence of a positive cosmological constant. We discuss the behavior of 
a quantized, massless and minimally coupled scalar field in the anisotropic stage. This scalar field 
is not the inflaton and hence does not contribute to the background dynamics. We focus on the 
quantization procedure and evolution in the pre-inflationary anisotropic background. Also, in this 
paper for simplicity the metric perturbations are not taken into account. The initial condition is 
set by the requirement that the scalar field is initially in an adiabatic state. Usually, in a quantum 
harmonic oscillator system, an adiabatic process implies the one where the potential changes slowly 
enough compared to its size, and the time evolution can be obtained from the zero-th order WKB 
approximation. In our case, such a vacuum state exists only for limited solutions of the anisotropic 
Universe, whose spacetime structure is regular in the initial times. In this paper, we call our 
adiabatic vacuum state the anisotropic vacuum. In the Kasner de Sitter model, for one branch of 
planar solutions there is an anisotropic vacuum unless fca 7^ 0, where fcs is the comoving momentum 
along the third direction, while in the other branch there is no anisotropic vacuum state. In the 
first branch, for the moderate modes, fca ~ k, where k is the total comoving momentum, the scalar 
power spectrum has an oscillatory behavior and its direction dependence is suppressed. For the 
planar modes, fca <Si k, in contrast, the direction dependence becomes more important, because of 
the amplification of the scalar amplitude during this interval of the violation of WKB approximation 
in the initial anisotropic stage. The qualitative behaviors in the Taub-NUT de Sitter models are 
very similar to the case of the first branch of the planar Kasner de Sitter model. 
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I. INTRODUCTION 



It is well-known that our Universe is filled with the thermal microwave radiation field, so-called the cosmic mi- 
crowave background (CMB). CMB follows the almost idealized Planck distribution of 2.73K. The intensity of CMB 
is distributed almost isotropically, but it contains small temperature fluctuations of order 10~^ to the background 
temperature, which contain much information about the history of the Universe. Recent measurements by WMAP 
satellite [ll-Q have shown that the observed map of CMB anisotropy is almost consistent with the Gaussian and 
statistically isotropic primordial fluctuations, which is nothing but the most important predictions of the inflation. 

However, after the release of WMAP data, several groups have reported that there seem to be a few anomalies in 
the CMB temperature map on large angular scales. The most well-known fact is that there seems to be suppression 
of the power of CMB fluctuations on angular scales more than sixty degrees 0]. There are other observational facts 
that imply the effect which induces the violation of the rotational invariance. More precisely, probably there are the 
planarity of lower multipole moments, the alignment between the quadrupole {i = 2) and the octopole (i! = 3), and 
the alignment of them with the equinox and the ecliptic plane 3 • There are other observational facts implying the 
large-scale anisotropy, i.e., odd correlations of ^ = 4 ~ 8 multipoles with ^ = 2, 3 multipoles 0, a very large, possibly 
non-Gaussian cold spot in 10 degree scale @, asymmetry of angular map measured in north and south hemispheres 
, even though some authors claim that there is no significant evidence for primordial isotropy breaking in five-year 
WMAP data [|. 

Indeed, to explain the origin of the anomalies, various solutions have been suggested, introducing a non-trivial 
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topology 0, a locally anisotropy based on the Bianchi type VII^ universe to explain the quadrupole/octopole planarity 
and alignment [loj . non- linear inhomogeneities 11 [ and assuming an elliptic Universe to explain the suppression of 
the quadrupole CMB power [l^l- More recently, in particular, models which introduce an explicit source to break 
the spatial isotropy, either during inflation or in the late time Universe, have been proposed, e.g., by the dynamics 
of an anisotropic energy- momentum component durin g in flation |13Hl(i | , by the large scale magnetic field [l7j , by the 
anisotropic cosmological constant [ll] or dark energy ]l9| . 

In this paper, we will revisit the possibility that such large scale anomalies may be the relics of the physics of 
preinflationary anisotropic Universe. Wald's no-hair theorem ensures that in the presence of a positive cosmological 
constant an initially anisotropic Universe exponentially approaches the de Sitter spacetime at the later time under 
the strong or dominant energy condition |20[. It implies that it is plausible that the Universe is highly anisotropic 
at the beginning. The cosmological perturbation theory in the pre-inflationary Kasner phase was formulated in Rcf. 
(21I [22! . They showed that in general in an expanding (planar) Kasner phase one of two polarizations of gravitational 
waves ^ is coupled with the scalar mode, but the other gravitational mode is decoupled. Of course, the tensor-scalar 
coupling vanishes in the isotropic limit. The latter mode can be amplified significantly before the onset of inflation. 
The instability of the gravitational wave mode ma y b e deeply connected to the unstable mode found by the Belinskii- 
Khalatnikov-Lifshitz (BKL) analysis [2^ (see e.g., 24| for review). In a contracting Kasner universe such an instability 
can be identified with the large scale anisotropic mode. 

The problem is how to set the initial conditions. In the standard infiation, the initial condition is set inside the 
Hubble horizon, where the effects of the cosmic expansion can be ignored and the adiabatic condition is satisfied. 
Therefore, to compare with the prediction from the standard inflation, it is natural to apply similar arguments to the 
case of the anisotropic Universe. Usually, in a quantum harmonic oscillator system, an adiabatic process implies the 
one where the potential changes slowly enough compared to its size, and the time evolution can be obtained from the 
zero-th order WKB approximation. In this paper, we follow this definition for the term adiabatic. In the standard 
inflationary models, an adiabatic vacuum is also defined in the same way. Note that an isentropic process implies 
(adiabatic) + (reversible) process and therefore is slightly different from our adiabatic one. We call our adiabatic 
vacuum state an anisotropic vacuum. In an expanding Kasner solution, there are two branches of solutions with 
the planar symmetry. Then, the answer to the question on the initial adiabaticity depends on the choice of the 
branch. In one branch, initially the expansion rate along the planar directions vanishes while that along the third 
axis is finite. In this branch, the initial spacetime structure can be seen as (a patch of) the Minkowski spacetime, 
and thus the corresponding anisotropic vacuum has the physical meaning. The tensor-scalar coupling mentioned 
above vanishes, and therefore the initial dynamics reduces to the three-independent harmonic oscillators. However, 
in the other branch the tensor-scalar coupling diverges in going back to the initial time. In this branch, there is no 
adiabatic state. Thus, we mainly focus on the first branch. For a given set of initial conditions, the power spectrum 
was investigated, resorting to the numerical ways in Ref. pll . [2^ . One of the important motivations of this paper 
is to aim more analytic understanding of the primordial power spectrum and in particular its direction-dependence. 
As the first step, we will consider a free scalar theory, which is the counterpart of the inflaton fluctuation. The other 
important purpose is to investigate the case of the other Bianchi model and to see how generic the prediction of the 
Bianchi I model is. In particular, we will consider the analytic solution of the Bianchi IX Universe, i.e., the Taub-NUT 
de Sitter solution. 

The paper is constructed as follows: In Sec. II, the background solutions in the Bianchi I and Bianchi IX models are 
introduced, say Kasner de Sitter solutions and Taub-NUT de Sitter solutions, respectively. In Sec. HI, we investigate 
the behavior of a massless and minimally coupled scalar field in the background of Kasner de Sitter solution with 
the planar symmetries. We discuss the existence of a well-defined anisotropic vacuum for each mode, in terms of 
the validity of the WKB approximation in the very early Universe and determine the initial mode functions. For the 
modes for which the anisotropic vacuum can be well-defined, we derive the final power spectrum and investigate its 
direction dependence. In Sec. IV, we repeat the similar discussions in the case of the Taub-NUT de Sitter solutions. 
In Sec. V, we close the article after giving a brief summary and discussion. 



Here, the tensor or scalar modes are defined in the isotropic limit, where perturbations are decomposed into modes on the maximally 
symmetric 3-space. 
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II. SET-UP 

We are interested in the cosmological solutions of the Einstein gravity with a positive cosmological constant in the 
absence of the matter field with the action 

s = iv/|i / d^xy/^dn - a) , (1) 



.2 

where g^^ is the spacetime metric, A is a (positive) cosmological constant and il/pi is the Planck mass. 

A. Kasner-de Sitter spacetime 

Assuming the spatial part to be flat and homogeneous, the solution to the Einstein equation is given by 

where 

a = In (ao sinh^/^ (^^^))' 
/3+ = |ln(tanh(^)), 

/3- = ^Vl-m2 1n(tanh(^)), (3) 

and |m| < 1 and H := -\/A/3. The solution initially behaves as a Kasner spacetime, followed by a de Sitter phase, as 
discussed in [HI, H^l- At the initial time r = 0, there is the Big Bang singularity. In the later discussion, we will set 
flo = 1, which can always be done by an appropriate rescaling of the spatial coordinates a;* {i = 1, 2, 3). 

In the later discussions, we will focus on the case with a planar symmetry m = +1 or m = — 1. In the branch of 
m = —1, in the early time limit r — >■ the scale factor in the planar direction e"+'^+ approaches a constant, while 
that in the x^ direction e"~^^+ linearly depends on r. Thus, in this branch, there is one dynamical direction in the 
early time. In the branch of m = +1, in the early time limit, e"~^^+ ~ r^/^ and e"~^^+ ~ t~^/^, and thus there are 
initially expanding planar directions and contracting x^ direction. Thus, both the directions are dynamical. 



B. Taub-NUT de Sitter spacetime 

A Bianchi IX universe has three Killing vectors satisfying = ekij^k- The corresponding one form = 

(— sin.T'^dcc^ + sinx^ cos x^ dx^ , cos x^ dx^ + sinx^ sin x^ dx^ , cos x^ dx^ + dx^) constitutes the spatial part of the metric. 
Therefore, the general spacetime metric of a Bianchi IX universe is given by 

ds^ = ~n{tfdt'+a{tf{x'f+b{tf{x'f+c{i)^{x'f- (4) 
For the case of a planar symmetry in the sense of Bianchi I model, a{t) — the metric Eq. (j4]) becomes 

ds^ = -n{tfdt^+a{t)'^{{dx^f+sm^x^{dx^f'^ + c{tf [dx^ + cos x^ dx^'^^ . (5) 
As a background solution, the Taub-NUT de Sitter solution is considered 

where 

A(t) -e + 2Mi + - ^ (^^ - " y ) ■ (7) 

In the case of A = 0, the solution reduces to the Taub-NUT solution. In the late time limit t — ?> cx), the spacetime 
approaches de Sitter (dS) with the Hubble expansion rate H = {klZfl'^. As for the Taub-NUT space, in the region 
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A{t) < where t is spacelike, the geometry is geodesic incomplete. In this work, we focus on the region where 
A{t) > 0. Requiring that there is no remaining anisotropy in the late time Universe, we choose = 3 and obtain 

A{t) ^ ^(^t* + 2e^t^ + 8Meh + e'^y (8) 

The function A(t) can be expressed as 

A(0 = ■^{t~to){t^ + tat^ + {2e^ + tl)t + itl + 2eho + &Me^)), (9) 

where to satisfies 

4 + 2eHl + SMfto + = 0. (10) 

For \M\ < there is no real root for Eq. (jlOp and for \M.\ > 2^/(3-\/3) there are two real roots (For 

M > 2£/(2a/3), to is negative and for M < -2i/{3\/3), to is positive). For \M\ < 2£/(3\/3), A is always positive, 
and the spacetime geometry approaches de Sitter in the past and future t — > ±00, while it is anisotropic during the 
intermediate time. In the later discussions, we will focus on the case of Ai > 2£/(3-\/3), where A(T) vanishes at two 
points, represented by t = to.+,to,-, where to.+ > to.- {to,+ = to,- for A4 ~ ±2£/3y/3). We set the null surface 
t ~ io,+ as the initial surface (In the later discussions, we omit the subscript "+"). Note that for to ^ 

By introducing the new coordinate T t — to, we obtain 

A(r) = ^{T^ + 4toT^ + {2e'^ + 6tl)T + 'i{tl + eHo + 2Me^)) 

= ^(T"^ + '^toT^ + + Gtl)T + 3tl + 2£\-^). (12) 

to 

In the new coordinate system, the spacetime metric can be written as 

ds^ = -n{TfdT^ + a{Tf ( {dx^f + sin^ {dx^f^ + c{Tf {dx^ + cos x^dx^^ ^ (13) 
To compare the time evolutions of two independent scale factors, it is useful to rewrite the metric 

ds2 = -n(r)2dr2 + e2"[e2^(^'((dxi)'+sin2a;i(dx2)') +e-4'^(^)(dx3 + cosxidx2)'] , (14) 



where a — (l/3)ln(a^c) and /? — (1/3) In (a/c). The dynamics of Bianchi IX Universe and in particular the Taub- 
NUT de Sitter spacetime can be understood in terms of an analogy with the classical mechanics, as is discussed in 
Appendix A. See [25| for the other Bianchi IX solutions which are coupled to a scalar field. 



III. SCALAR FIELD IN THE BIANCHI I MODEL 



We consider a massless, minimally coupled scalar field propagating on the background spacetime discussed in the 
previous section with action 




Note that this scalar field is not the inflaton and hence does not contribute to the background dynamics. We focus on 
the quantization procedure and the evolution in the anisotropic stage. The effects of a scalar field onto the geometry 
of the anisotropic Universe were discussed in [2^ . 
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A. Equation of motion 

Explicitly, the scalar action is written as 



2 

Introducing the Fourier decomposition for the mode function 



i f drdPxe'" (0^^ - e.-'^'^h'^ d^<j)dj(l^ . (16) 



where 7 = (fci,fc2,fc3) denotes the mutually orthogonal components of the comoving momentum, the action can be 
reduced to 



1 



V drj 



2 



S^ = -j d'kdr, ( ^ ) - io'^x' , (18) 



where Xi = e'^ff'-y drj ^ e "'dr, which corresponds to the conformal time in the isotropic case. The frequency 
squared is given by 

where H — \/ ^ is the Hubble parameter for a large r and the momentum dependent part is 



1 + Vl + ) ^ ( 1 + \/l + e^" ) ' ^ g4(m+l)a 



f]^(t) = ^ ' k\ + ^ ' fc^ + kj. (20) 

g2(m-2+V3(l-™^))a g2(m-2- ^3(1-™^))^ (l + V l + e^") 

The equation of motion for each Fourier mode becomes 



d Xi 
drf 



+ <X7 = 0. (21) 



In the following discussions, we will focus on the background solutions which have the planar symmetry, m = —1 
or m = +1. For the branch m = — 1, in the limit a — > 0, 

^ +(-^^+372^)- 
where A:^ := fc^ + fc| + /s| — k^ + k^. In the later time limit e" — > oo, 

= -2^^^^" + - + 1^ + 0(e--) , (23) 

which approaches the one in the ordinary de Sitter universe for large scale factor e". The case of the branch of 771 = +1 
will be discussed in subsection II. G. 

B. Quantization and power spectrum 

The canonical quantization of the scalar field is done in the standard manner: 

d^k(ui,aw + Kaiy (24) 

where [cki , a^^] ~ '^(''1 ~ (others are zero) and Uk = e*''^^k/ {^tt)'^^^, and the mode function in the time-direction 
can be normalized as ^kC^t'/'k ~ {dt4''k) 4't. ~ V^'^"- The power spectrum is defined by 

(010^10) I d\nk I ^P,, = (25) 
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Note that in contrast to the case of the standard inflation the direction dependence would be included in the power 
spectrum. The vacuum is chosen at the initial anisotropic era: r — >■ 0—. Usually, in a quantum harmonic oscillator 
system, an adiabatic process implies the one where the potential changes slowly enough compared to its size, and the 
time evolution can be obtained from the zcro-th order WKB approximation. In this paper , we follow this definition 
for the term adiabatic. In the standard inflationary models, an adiabatic vacuum is also defined in the same way. 
Note that an isentropic process implies (adiabatic) + (reversible) process and therefore is slightly different from our 
adiabatic one. In the anisotropic spacetime an adiabatic vacuum state can be found only in the special solutions of 
the anisotropic Universe, which are regular in the initial times. We call our adiabatic vacuum an anisotropic vacuum. 



C. WKB solution for m = -1 



From now on, we will focus on the branch of m = — 1 except in the subsection II F. By changing the coordinate 
dt = dr/e^", which leads to 

e" = sinh^/^(3i/T) 



sinh^/^(-3iJ0 ' 

one obtains the equation of motion in the form of a time-dependent oscillator 

When r changes from 0+ to +00, t changes from —00 to 0—. The frequency squared becomes 



(26) 



(27) 



where r± := k±/k and x{t) = 1 — e^^* = e^®"(-\/e^" + 1 — 1) monotonically varies from one to zero as time t increases 
from negative infinity to zero. 
The WKB solution is given by 



1 



»WKB — 



2n(t) 



exp 



- i / dt'n(t') + ill) 



(28) 



where '0 is a phase factor and i/|io| 3> 1 and 17 satisfies the nonlinear equation 



2_ 3^ 



5i(t)2 = n(t) 



The WKB wavefunction is valid up to the order of the correction term if 



e(t) 



dOf{t) 
dt 



n^(t) 



H 



1 - x(t) 



^ (x/2)y^{l~rlx(t)Y'^ Vl-^i-^W 



1 < 1, 



(29) 



(30) 



where e play the role of the adiabatic parameter and is plotted in Fig. [T] for several different values of r^. At a; = 1, all 
the curves converge to zero and therefore WKB approximation is always valid at the beginning. Thus, it is possible 
to define an anisotropic vacuum. In the case of m = —1, the initial spacetime metric can be approximately written as 



ds^ 



-dr' +d{xy + d(xy +T'd(xy, 



(31) 



where unimportant constants are absorbed by the rescaling of . It is straightforward to sec that this metric represents 
(a patch of) the Minkowski spacetime: 



ds' = -dU' + d(x')' + d(x')' + dR 



(32) 



by the coordinate transformation R = rsinlix^ and U — rcoshcc^. Thus, it is possible to find the positive frequency 
modes and the associated vacuum state, for an observer moving toward the increasing r-direction along a constant 
x^ curve on the (t/, i?)-plane. This corresponds to our choice of the vacuum at the beginning time. 
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FIG. 1: The behavior of e(t) with respect to x. Here, = 0.999,0.99, .9, .5, and .1 respectively from the top and we choose 
H/k = 0.05. The WKB approximation is valid if e{t) < 1. The matching time t, is almost independent of k. 



Near x ^ 1, the condition is always valid for all momentum modes except for the mode with r± = I. The condition 
is also valid for high momentum modes if the scale factor e" satisfies e" <C -g^ i.e., —Ht ^ ("f")"^- On the other hand, 
the large e" approximation (|23p is valid if e° 3> 1. For high frequency modes, if we choose an intermediate time 
satisfying 



1 « e"(**) « 4 
H 



(33) 



and choose the WKB wavefunction for t < and the large e" approximation, for t > the approximation will be 
good and in fact its accuracy is almost independent of the choice of the intermediate time as long as it satisfies 
Eq. Therefore, we choose 



(34) 



and the adiabatic parameter at t = t,, is order of 0[H^^^ /k^^^) . For modes with k/H ^ 1 or k < H, the 
approximation is not valid. The WKB wavefunction is expanded up to an enough adiabatic order, to validate our 
matching scheme with the solutions in the quasi-de Sitter region. The assumption of fc ^ i7 is justified if we assume a 
slow rolling inflation for our universe. Note that however, as we will discuss later, for a given class of inflation models, 
the modes oi k > H arc observable in some cases. 

Near a; ~ 0, the WKB solution at the zeroth adiabatic order takes the form: 



OWKB 



exp 



(35) 



Note that the next order correction at e"'-**-' = \/k/H is 0{H/k) since :— x{t^. 
The value of the frequency fi* := ilfc(i*) is 



2(i7/A;)3/2(l_ (ii-/fc)3/2 



H 



3/2 



The integral on the exponent is evaluated as 

ft 



k I dt 

J to 



2^/Vl-ri(l-e6g*) 

(1 — g6/ft~)2/3 



L — X xP-/^ 



(36) 



(37) 



where e := e^^*" is some cut-off scale at an initial time. The integration can be executed explicitly to give the Appell 



hypergeometric function of two variables, 

dx (l-rix)i/2 



1 



r;2/3 



2.f(l-ri.x)3/^ 3 2 5 ^ _ ^ l_ 
3(1 -r2^) ^4' 3' '2' ^ ' 1 
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D. Matching 

In the (quasi-) de Sitter region, (— t) ^ H^^ , the equation of motion (p6)) becomes 



(— 



)0 = O, 



whose solution becomes 

where the positive and negative frequency modes are given by 

ik. 



(38) 



(39) 



This approximation is vahd for e" 3> 1. The WKB solution is now matched to the one in the quasi-de Sitter region, 
at i = i*. Now, we should determine the adiabatic order where the matching is done. As an example, we post the 

frequency n at U by using = sinh"^(e-3"(**)) = sinh~\(i7//s)^^^) and Eq. 



H 



l + 2(--ri) 



3/2 



(40) 



As seen here, the first direction dependence appears at the third order in e'^ and wc need the WKB approximation up 
to the fourth adiabatic order. Then, the WKB frequency is given by 



The matching condition is given by 



A+(^+(t^) + A_0_(t,) = 1 + 



t. 



2k J JWZ 



A+(j)+'{U) + A^(j)J{Q = -i 
Now the solution is exact up to the adiabatic order of e^. The phase factor is given by 



H, H f H 

^'-2k^-lTV-2k 



<f>* ~ exp 



I I fldt + iip 



It is straightforward to confirm that the coefficients A^ and A^ satisfy the normalization condition 



U4 



Ml 

2fc3 



As the reference, in the standard isotropic inflation, the normalized mode function is given by 



V2fc3 \ -tl / 



1/3 



(41) 



(42) 



(43) 



(44) 



(45) 



E. Power spectrum 

In the limit of the later times t 0—, the wavcfunction behaves as 

4>fin = -(A++A_). 
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Explicitly, the wavefunction in the later time to the order becomes 








^)( 


J) 






(1 


-rl 




Then, the power spectrum Eq. (|25l) in terms of the adiabatic order 0{{H/k)^/^) is given by 




3/2 



2H\ k fH 

- cos2^- + - 



(46) 



As a reference, we present the power spectrum from the standard inflation obtained from Eq. (P5|) and (|35 



p(0) ^ 



2 

47r2- 



There are two modifications: One is an oscillatory behavior of the power spectrum, which seems to be observed in 
the CMB angular power spectrum around the first acoustic peak ^ = 10 '--^ 100, and the other is an overall change of 
the amplitude given by (1 — H/k). In the above expression by taking k ^ H, the power spectrum seems to vanish. 
This would be consistent with the suppression of the angular power spectrum of CMB fluctuations on large scales if 
the anisotropy of the present Hubble horizon scale corresponds to the wavenumber of the order k > H. For larger k, 
it approaches the scale invariant form with a small oscillatory behavior. The direction dependence appears only on 
the oscillatory term. 

As discussed above, our WKB approximation is valid only for modes with k/H ^ 1. Indeed, one may ignore 
the contributions of modes such as k ^ H or k < H. One may not observe modes whose wavelengthes are longer 
than the present Hubble scale. Thus, there is the critical wavenumber given by fccrit = o,oHo, where H^^^ ^ lO^^cm 
is present day's Hubble scale: one observes the modes of A: > fccrit- If it is possible to take fccrit > H, our WKB 
approximation is justified for all the observable modes. On the other hand, Qq = (T^/Tojaf, where Tr, Tq sa 2.7K 
and Of are the reheating temperature, the CMB temperature at present and the value of scale factor at the end of 
inflation, respectively. One finds the relation 



7 TT ^jT-R jj N Th Ho tj 

fccrit = aoiJo = -^7^ Ho = e"——H ■ 
J-o J^o n 



AH, A := 



Tr Vbe 

^ To ^1/2' 
mf 



(47) 



where N = \n{a[/a[) ~ In at is the number of e-folding of inflation, Vint and Vde ~ (10 ^eV)** are the values of 
inflation and dark energy potentials. Note that H w ^ {8'K/3M^)Vinf is the Hubble constant during inflation and 

Hq ~ w (87r/3Afp)VDE is that of the present Universe. One can estimate the value of A as 



A 



Tr 



lOi^GeV 



10i4GeV 



V 



1/4 



inf 



Thus, for inflation with GUT or a lower energy scale and if e-folding is around N > 64, it may be possible to realize 
A> 1, leading to fccrit > H. On the other hand, for sufficiently long duration of inflation, say N ^ 64, the observable 
modes are only those with k ^ H since A 3> 1, and therefore it seems to be very hard to find any deviation from the 
scale-invariant spectrum. In other words, to find the information about the preinflationary anisotropy, the number of 
e-folding must be tuned. 



F. Planar modes, r± < 1 

As seen in the case of = 0.999 in Fig. [U namely for modes on a plane, there appears a region where the WKB 
approximation may not be valid near x ^ I. We divide the time into three separate regions divided by the times ti 
and i*. In the region ti < t < the WKB approximation is valid. For other two regions, we may find an appropriate 
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approximate solution. Note that in the case of rj^ = 1 exactly, the adiabaticity parameter diverges in the limit of 
i — ^ — oo and there is no anisotropic vacuum state. 
The frequency squared for x ~ 1 becomes 

n'- = 2^'\kl + kle^"') + 0{ee^^"'), for Ht « -1. 

The corresponding solution for Eq. \26\ is 



QH sinh(7rg3) 



J_,,3((7ie^^*), foriJi<-l, 



(48) 



where 



2^/3 1^31 

93 = ^rr ^ 9J 



3H ' 'iH 
Here we chose the solution so that it becomes an incoming wave form: 

1 



V2-22/3|fc3| 



exp 



I - 22/3i|/c3|< + fort 



— oo, 



where the initial phase becomes e*^ = [r(l + iq3)/r(l — 153)]^^^ (^)*'^- We choose the momentum dependent cutoff 
time ^1 so that both of the solution (l48l) and the WKB solution are equally valid. The WKB solution is valid if e < 1. 
For a:; ^ 1, we check the validity of the WKB approximation, 



3 • 21/3/7- e^Ht 



3 • 1 



k (r2+e6^*)3/2 k 

where in the second equality we assume e3^* ^ since we want to be placed over the first peak in t as seen in 

Fig. [1] For the WKB approximation to be valid, we should choose the time so that e'^^*^ > '^'^ ^ ^ . For example, we 
may choose the time ti satisfying 



„3ffti 



22/3 V fc' 



(49) 



where both of the WKB approximation and the solution (|48|) are equally good. Note that the time ti is direction 
independent. At time ti, the assumption e3^* ^ r3 constrains the orthogonal momentum value to 



(50) 



For planar high-momentum modes, the argument of the Bessel function, gj^e^^*^ 
Therefore at ti, we may employ the asymptotic form for the Bessel function: 



lim Jia{x) 



During ti < t < t^,, we use the WKB approximation where is given in Eq. ([M)) . The WKB solution is 



-g-, is very large. 



(51) 



V2O 



exp • 



{-il dt'n{t')] + ^=cxp\i I dt'n{t')] 

^ Ju ' ^/2^l Ju ' 



(52) 



For t ^ t^, the WKB approximation is not valid. However, for t > t^, the asymptotic solution (|38)) is valid. 
At ^1, to the zeroth order in H/k, the matching condition becomes 



Mti) + i 
01 (ti) - i 



^2sinh(7rg3) 

—i 

^2sinh(7rg3) 



H^«2\ -1/4 



cos V' — M 1 + 



sin'0 



f7«2\ -1/4 



cos + « 1 + 



2 \ 1/4 



sin'0 
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where n 



2,2 , 2^/-^ ,2 



the coeffLcient B± becomes of the zeroth order, 



(i + f¥) 



and ip 



2 '4 



/A; if 



Since g| < k\/{kH) as in Eq. ([50]) . 



^4 



X e 



vr 



\/l - e-^^r^s 

For < t < 0, the isotropic solution is valid and 

03 = C+0+(O+C_ 



-27rg3 



X e 



(53) 



(54) 



The primordial power spectrum is proportional to the square of the size of C-f + C_ . Matching this solution with the 
WKB solution ([5^ at the boundary , we have 



C- = 



2,fc3 [02(t*)(0-(t,) - Mt*)) - (l>'2{t*){4>-{U) - 4>+{U))\ 

Using Eqs. and we get, in the limit g| < k\/(kH), 



C+ + C- 



H 



{-(1 + gin( J - 4>(fc)) + z(l - e-'^'?--^) cos( J - 4>(fc))} , (55) 



where $(fc) = y |^ + J^*^* dt'fl{t') + y ji and we ignore terms of 0{^jH/k). 
The primordial power spectrum to the zeroth order is given by 



(56) 



1 + icosech^^ (1 - sin(2$(fc))) 



If one reduces keeping k constant, the direction dependence of the power spectrum becomes un-negligible at 



93 



2/7r or ^3 



3.25/8 



H <^ k. Note that this result is quite different from the general expectation [ij] that the 



direction dependent term would be (approximately) scale invariant and second order in fca. The phase factor $(fc) is 
dominated by J^* dt'Q{t') and is becomes 



2*(fc) = 2\/^ + 2 




dt'n{t') 



22^VHXl/3)_fc 
3r(5/6) H' 



G. The case of m — +1 



In this subsection, we briefly discuss the case of the other solution with a planar symmetry, m ~ +1. We will see 
that the WKB approximation can not be applied at the beginning and there is no anisotropic vacuum. In this case, 
it is appropriate to change the variables as (/) = a~^^'^^ and dt = a~^du, to rewrite the equation of motion into the 
form of a time-dependent harmonic oscillator: 



where 



du'^ 



+ E{uf $(u) = 0, 



(57) 



n2 " 



-2o2 



4a2 



Za 



,ri(l + y/l + a{u)^) + ajufjl - rj) 
(1 + v/l + a(u)S)4/3 



(58) 



and the prime denotes the derivative with respect to u. The coordinate u is defined in the range < u < uq 
•\/7rr(l/3)/(37Jr(5/6)) « 1.4027J~^. In order to check the validity of the WKB approximation, one needs another 
coordinate transformation such that Hv = ln(it/ito) and $ = e'"^^^'^, where v runs from — oo to zero, and then the 
equation of motion becomes 



-dv"^ 



* = 0, 



(59) 
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where 

■.= H'uy-E{u{v)r-^. (60) 

The validity of WKB approximation is determined by the parameter €„ :— [(d^^/du)/^"^ | . For t — > — oo, u — >■ and 
V — > — oo, then a w Hu, and one can evaluate « k^r^H^UQC^^'" and e„ « 2e~^^ /k±^UQ . Therefore, the adiabaticity 
parameter exceeds unity for v <C — (1/i?) In (/c^uq). In such a case, one cannot find a well-defined adiabatic vacuum 
state. 



IV. SCALAR FIELD IN THE BIANCHI IX UNIVERSE 



In this section, we investigate the primordial spectrum of a massless and minimally coupled scalar field in the 
Bianchi IX Universe. Such an investigation will give a deep insight about how generic the results obtained in the 
case of the Bianchi I model are. One reason is that recent WMAP data seem to favor a spatially closed Universe ^, 
and the Bianchi IX Universe exactly has that geometry. The other important reason is that, similarly to the case of 
the Bianchi I model, in the Bianchi IX model there is an exact analytic solution, the Taub-NUT de Sitter spacetime 
whose metric is given by Eq. ([5|). The spacetime geometry is a Taub-NUT spacetime in the initially anisotropic era, 
and approaches de Sitter solution in the later time. Our discussion will come along the same line in the case of the 
Bianchi I. Note that as discussed in Sec. II, we restrict the value of to be > 21/ (3-^3), for which A(f) vanishes 
at i = <o < (r = 0). For such a choice, one can naturally choose the initial time to be T = 0. 



A. Scalar field on the Bianchi IX Universe 



The equation of motion of a massless and minimally coupled scalar field, 00 = 0, explicitly reduces to 



1 



1 



(-l + e«^)(-Ll) 



= 0, 



(61) 



where the angular momentum operators are defined by 



1 4^r.ntT.^ 



+ cot X 



1 



dx^ ' sin^a;! 



\dx'^J \dx^J dx'^o 



The general classical solution can be decomposed by the basis of harmonic functions (j}{T^x^) = 
(/3^(T)y^(x^, a;^, x"^), where 7 = {J,K,A1) represents a short-hand notation of a set of quantum numbers 
associated with the spatial directions (27| . The harmonic function on the three-space can be constructed by 



0-y(a;^), where Qjkm satisfies the ordinary differential equation 



cot X 



+ ~ 2 cos x^ KM 



.d{x^f d{x^) 
The harmonic function is normalized as 



• 2 1 
sm a;^ 



dx^dx^dx^ sin x^Yj{x^ , x'^ , x^)Y* {x^ , x'^ , x^) = ^7,7'. 



The equation in the time-like direction is given by 

d^ 



= 



J( J + 1) - K' 



n'dTi''"' ' ' e2("-2« 
For a given vacuum |0), the quantization of the scalar field can be done in the standard way 



(62) 



(63) 



^ Simultaneous constraints on the spatial curvature and the (constant) equation of state parameter of dark energy by the five-year WMAP 
data seem to favor a positive curvature, although the flat Universe is also consistent at 95 percent confidence level 
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where the creation and annihilation operators, ajj, and a^, have the properties a^|0) = and 



7t 



0, and satisfy 



the commutation relations [a-yjjajj,^] — (^71,727 and [a-y^jO^^] = [aj^j^jflt^] = 0. In our later discussions, a vacuum |0) 
will be chosen in order for the initial modes to be adiabatic. The mode functions satisfy the Wronskian normalization 
condition 



For the case of the discrete spectrum, the power spectrum can be defined 



(O|02|O)=^P^, := 



7 

It is convenient to define the new time coordinate as 



2{T^ + 2ioT + P0 + I) 



2£A T{T^ + 4X2 + (6*2 + 2)r + 3*g + 2to - l/^o 



(64) 



(65) 



(66) 



and variable a.s ip = e " ^Xi where T := T/£ and to to/£ arc dimcnsionless. Then, the equation of motion is 
rewritten as 



da; 2 



+ r!(.T)2 x{x) = 0, 



where the frequency squared is given by 

n^{x) ■.^nl{x) + (^j{j + i)-K^y-''^ + 

Here the momentum independent part is 



(67) 



(68) 



n2(T) 



/da d/3\2 /d2a ^2^. g2a-4;3 ^2^, ^2^ 2da d/3 Idnwda 



dx^ J 



\dT^ dT^J \dT dT 



.dTJ \dT dT) 



TF{T) + 4toT2 + 6P0T + 2T + 3if, + 2io - l/k 



F{T)AiT) 



4 (r2 + 2ioT + + 1) 

where A is defined in Eq. ([7]) and 



i^T^ + 2ioT + il + iy 



(69) 



F{T) 



2T'' + I2TH0 + 2,QT% + 5r* + 4ir^tg + 22TH0 

IhilT + 2AtlT +7ioT-2Z + 3il + + P0-I. 

to 



to 



■2-7-2 



(70) 



For a small T, 



and for a large T, 



ii-3Po? ^ 
4fo(l + il) 



(-2 + 25t^ + 6ig - 66i^ - 18ig + 27?g) 
4tg(l + ig)2 



For a large T a;(T) = — 1/T + 0{T^), and for a small T the relation becomes 

1 - 34 To 

where To {M / i) can be determined if we integrate out the equation (|66| numerically. Note that 

_ 4A/^2 ^ r(r^ + 4for2 + (6*1 + 2)r + 3tg + 2io - i/ip) 



(71) 



(72) 



ir' + 2Tto + Po + iy 



(r2 + 2rio + + 1)2 
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B. WKB solutions, matching and scalar power spectrum 



We discuss the primordial spectrum in the Bianchi IX Universe. The new coordinate x runs from — oo to 0—, where 
T — > 0+ and T — >■ oo correspond to a; = — oo and x ~ 0—, respectively. The high angular-momentum limit (J — > oo) 
corresponds to the high momentum limit fc 3> in the Bianchi I model, since J is dimensionless. 

a. Adiabatic parameter e and the WKB approximation: We define the adiabaticity parameter, 

_ 1 dn^ _ e"-^^ 1 dn'^ _ 2£A 1 dn^ 

eK^o--\e\; ^-if^p^^.- „ \n\3 dT ^ + 2aT + + \n\3 dT ' ^ ' 



0.3 



0.2 




FIG. 2: The behavior of e(T). The horizontal axis represents T. Here we take K = 2 and J — 30. The parameter to is chosen 
to be to = —0.5, —0.3, —0.1, —0.05, —0.03, —0.01 respectively from the top. As seen in this figure, the large J approximation 
gives the best result at T ~ Vl and gives a good approximation even for time T ~ J/2. 



We focus on the case of a small T (|T| ^ 1). The value of e vanishes at T = and varies with time as 
_ [4(,/(J+ 1) - + 1] / ^,~2,W 3[4(J(J+1) +l]2_^32A' 



(l + 0(?„)) T -^-^ ' ~ >^ ^ (l + 0{tl))r- + 0{T'). (74) 



8(-fo) ^ ^ 64(-to)3if5 

As in the m — —1 branch of the planar Kasner-de Sitter model, the WKB approximation is always valid at the 
beginning and therefore one can find a well-defined anisotropic vacuum state. This can be understood as follows. 
From Eq. (13), with use of Eqs. (6) and (12), it turns out that in the T — > 0, the initial metric functions behaves 
as n oc T^^/^, a oc constant and c oc T^/^. After defining the proper time dr ~ n{T)dT « dT jT^I"^^ one finds 
a oc constant and c oc r, which is a patch of a static spacetime as Eq. pip . 
As T increases e(T) has a maximum value 

[4(j(J+l)-/r2) + l]' 1 



4(-io)/^ f3[4( J( J + 1) - i^2) + 1] 2 _^ 32^2 



12if(-to) 



where the last step is for large J at 

4A-2(-to)[4(j(,/+l)-/v2) + l] 
3[4(J(J + 1)- A'2) +l]V32/f2 ' 

and then start to decrease toward negative values crossing zero. For modes satisfying emax < 1, we may use the WKB 
approximation around T = Tmax- Therefore, assuming J ^ 1. the WKB approximation is valid at the initial period 
of time for modes satisfying. 



(75) 
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The subsequent behavior of the adiabatic parameter can be described by the large J approximation. 



2to^j{j + i)T{T^ + 2kr + ii + 1) + 4fgr2 + eiir + 2r + afg + 2i„ ~ i /to 



Starting from a positive value around T ~ 0, the explicit value of e decreases and vanishes at T = 1/^/3 — to- In small 

to limit, it continually decreases to a negative number of the order —1/ fp) and then bounces back to zero as T 
goes to infinity. Therefore, the WKB approximation in this intermediate region is valid if 



This large J approximation, however, is not satisfied both for small and large T limits. As for a small T, Eq. ([71)1 
is a good approximation. For a large T, the highest power terms in T will become comparable to the large J 
contributions, and the adiabatic parameter can be written as follows: 

V2[l + 3(-to)-i(l + Po?{JiJ + 1) - K^)/T'] 



|2J(J + l)/(r2) - 1 



1 3/2 



The numerator will be dominated by the J dependent term for a time T smaller than Tm = {J{J + l)/io)^^^- On the 
other hand, if T > %n, we may neglect the J dependent term compared to unity. The denominator vanishes at time 
Tc = v/2J(J+ 1), leading to the divergence of the adiabaticity parameter. For T < Tc ''^ J, the first term rules and 
for T > Tc, the constant (—1) rules the denominator. Noting these behaviors, the order of e will be minimized at time 
7" ~ Trm where the accuracy of the WKB approximation will be maximized. At a time T satisfying Tm < T Tc, 
the adiabaticity parameter can be approximated to be 

7^ 

(77) 



2[J(J + 1)] 



3/2 ' 



which is much smaller than unity. 

The adiabaticity parameter then goes to negative infinity at T ~ Tc (for a large J), and then bounces back to 
a finite negative value —V2 as T goes to oo. As T — > oo, e approaches —V2 which is larger than one and the 
WKB approximation is broken there. In this region, we use the large T approximation. It is good to solve the later 
time evolution equation (j62p in T coordinate rather than in x coordinate since the order of the frequency part is 
transparent: 



J{J + 1) - Kl ^ ^ 4J( J + 1) 4?o , 10^1 



=2/3 



1 - ^ + — ^ + OiT-') . (78) 



7-4 V T T^ 



The first T term gives the asymptotic solution x± in Eq. ([5T|) below and the second T ^ term gives corrections. 
Therefore, at a time T ^ 1(^ to)- the error of the solution x± is of order 

T • 

The matching time T* is determined by the condition that this accuracy is the same as Eq. ([77]): 

M-io) ^ T' 

T 2{J{J+I)f^' 

We get the matching time 

%^2[JiJ + l)fi-io)^/\ (79) 

In summary, the WKB approximation is valid for T ^Tc for modes satisfying Eqs. ([75)1 and ((75)) . However, it fails 
to be satisfied for T > Tc and we use the large T approximation. These two solutions are matched at time T, and 
the accuracy of the whole solution is of 0{{—io)'^^^J~"^^^). 
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b. Moderate modes: If |to| is smaller than 0(1) {Ai is not so large), most high momentum modes except the 
modes K < \, satisfy Eqs. ([75]) and The main contribution to the power spectrum comes from these modes. 

For them, we use the WKB solution, which is given by 

XWKB(a;) = J— exp I - z / + #|, t < tc- (80) 

y/2\l{x) >- JxQ ' 

Approximating CP' with 4J(J + 1)/T'*, the later time solution is given by 

X{x) A+x+{x) + A^x-{^), 

1 -i(,/(,/+i))' 



V/2(J(J + 1))'/' V {J{J + l)Y''x, 



X-(x) := ^ e<"(-^+^))' -1 , (81) 

V/2(J(J+1))^/^ V (J(J+l))^/^xy 

where x+.xX- ~ X-,xX+ = ~* normalization condition implies the relation |j4+P — = 1. 

The matching at x = a;* determines the coefficients 

A+ = -^^{x'-{x*)+^^*X-ix*)), = ^^{x+ix*) + ^^*^+(^*))' (82) 

where <I>* is an unimportant phase factor in WKB solution. The power spectrum is given by 

1 



JKM 



2£2(J(J+1)) 
1 



3/2 



A+~A. 



The later times the frequency squared at the present order becomes 



v/J(J + l):r,c.- [^[^t^,W4 ' ^■/(■/+l)l/l- , -^v/^(J+l). 

Therefore the power spectrum, including the first non- vanishing direction dependent correction, becomes 



^ , 8{-i,f\m{[JiJ + l)f\-io)-^/^y 



2^2(J(J+l))3/2 ^ [J(J+1)]'^' 



(84) 



As the reference, it is useful to show the scalar power spectrum in the standard de Sitter invariant vacuum is given 
by 

^j"^'' ^ 2£2(J(J+l))3/2 (85) 

and thus multiplied by the factor ^ J^, which is proportional to the density of state for a given J, the spectrum 
becomes scale-invariant. Thus, in the resultant power spectrum for the moderate modes, the deviation from the 
standard spectrum is suppressed by the factor 0{J^'^/^). 

c. Planar modes with K < - . ; In the case of planar modes, the WKB approximation fails to be satisfied at 

io } 

^^2/ ^ \ .... . . . ~ 

T ^ Tmax ^ — 37^- Around these initial times, we try to find an approximate solution for the times t to. We 
need to expand the frequency squared up to first order in toT: 

(K^+ '^'^'}}'-''!'f ], (86) 
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where we ignore 1/io and K relative to J. Noting '^'■'^"^jf^ "-^ ~ y , we get the solution to the differential equation (|62p . 



y (1 - 3ig) smh (TrgfKj 



where 



l-SlJ ' V(l-3'l)(l + 'l) 

and we choose the normalization of the solution so that it becomes an incoming wave in x coordinate at x — >■ — oo, 
noting the relation ([7^ between T and x. If we use the solution (|57| up to a time 71, the accuracy of the above 
solution will be of order 7i(— ^o)- 

On the other hand, the WKB solution is still valid at a time 71 satisfying -^^ — ^ 7i ^ 1. At this time we 

may also use the high J limit and the adiabaticity parameter becomes e ~ 2(-t yl'^JT^i'^ ' ^^'^ matching time 7i 
between the initial exact solution and the WKB solution can be determined by setting the accuracies to be the same: 



Therefore, the matching time becomes 71 = — ^ = t ■ The WKB solution becomes 

^ 22/3(_to)J2/3 



We use this solution during the period of time Ti < T < %■ The coefficient B± is determined to be 

= ^ ^'"^ L m) ± - ^^^^^ ^ 1 + ^ ,..(A.rV--,) (90) 

^ 2x/2 l^-^i^^i^* 2(-fo)f^i rfr y 2 V l-e-2-9K ' ^^"^ 

where we write down the zeroth order term only. 

After T > T*, we may once again can use the asymptotic solution x± to compare it with the WKB solution xwkb • 
The late time solution is given by 

X{x) = c+x+(a;) + c^x-{x). 

At r*, fi* = J(J + 1) and (J( J + l))i/2x, ~ -J^'^/{-iof'^. 

The primordial power spectrum is proportional to the square of the size of c+ — c_ . Matching at the boundary , 
we have 



dx- dxwKB \ . ( dx+ dx 



c+ = -M Xwkb—, j X~ , c_ = z xwkb 



iWKB 



dx dx J V dx dx 

If we keep only the zeroth order, the power spectrum becomes 



JKM 



c_2 = °- — ^ 5 5 cos2$ , 91 



2^2(j(j+i))3/2' + 2£2(j(J+l))3/Hl-e-2-9^^- l-e-2-«K 



where $ = + l))^^^^;* + /J^' fidx + AjTi^^ — f and = A{~-io)K/ (1 — Stg) contains the direction dependence. 
Thus, the explicit direction dependence is not suppressed. Note that ^ J rules the phase factor $ for larger J. The 
K = limit is not well defined, however, this mode is not in our concern since the wavelength of such a perturbation 
mode is beyond our Hubble horizon. 
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V. CONCLUSION AND DISCUSSIONS 

111 this article, we considered the quantization of a massless and minimany coupled scalar field in the Universe, which 
is initially anisotropic and approaches the de Sitter spacetime. The motivation to consider the initially anisotropic 
Universe is two-fold: The first motivation is that even if the current Universe is almost isotropic, it does not mean 
that the Universe is isotropic from the beginning. In fact, Wald's no-hair theorem ensures that in the presence of 
a positive cosmological constant an initially anisotropic Universe exponentially approaches a de Sitter spacetime at 
the later time under the strong or dominant energy condition. Therefore, it would be more generic that the initial 
Universe is anisotropic. The second motivation comes from the recent observations by WMAP satellite. WMAP 
measured the temperature fluctuations of cosmic microwave background (CMB) and almost confirmed the predictions 
from the inflation, during which the Gaussian and statistically isotropic primordial fluctuations are produced. But 
after the WMAP data were released, several groups have reported the so-called low-^ anomalies in large angular power 
of CMB fluctuations, e.g., the suppression of the power of quadrupolc, the planarity of quadrupole and octopole CMB 
maps, the alignment of the preferred directions of quadrupole and octopole moments, and so on. They may not be 
satisfactorily explained by the standard isotropic initial state, but may be done by the direction-dependent primordial 
fluctuations. 

In this paper, we considered the gravitational theory composed of the Einstein-Hilbert term and a positive cosmo- 
logical constant. For simplicity, we assume that the late-time inflationary stage is exactly described by the de Sitter 
solution and ignored the dynamics of an inflaton field for simplicity. We considered two kinds of initially anisotropic 
Universe, say, the Bianchi I and Bianchi IX models. In each model, there is an exact solution. In the Bianchi I 
model, we considered the Kasner-de Sitter solution, in which the spacetime geometry is initially a Kasner spacetime 
and approaches the de Sitter spacetime. In the Bianchi IX model, similarly we considered the Taub-NUT de Sitter 
solution, in which the spacetime geometry is initially a Taub-NUT spacetime. Note that in the case of Bianchi I, we 
focused on the case that the Universe has an exact planar symmetry, which is isotropic along two of three spatial 
axes. In the Bianchi IX model, the exact solution has only two independent scale factors. 

After giving the background geometry, we investigated the spectrum of a scalar field, which is the counterpart of 
the inflaton fluctuation. We discussed how a massless scalar field is quantized in the initial anisotropic stage. The 
qualitative behaviors of the resultant spectrum in both the Bianchi I and IX models are very similar and therefore, here 
we summarize our result mainly focusing on the Bianchi I model. We found that in the case with a planar symmetry, 
there is a well-defined adiabatic vacuum unless fcs 7^ 0, where fca is the comoving momentum along the preferred 
direction. As we mentioned previously, in a quantum harmonic oscillator system, an adiabatic process usually implies 
the one where the potential changes slowly enough compared to its size, and the time evolution can be obtained 
from the zero-th order WKB approximation. We followed this definition for the term adiabatic. In the standard 
inflationary models, an adiabatic vacuum is also defined in the same way. In our case an adiabatic vacuum state, 
called an anisotropic vacuum in this paper, was found only in the special solutions of anisotropic Universe, which 
are regular in the initial times. It was shown that for the moderate modes, ~ fc, where k is the total comoving 
momentum, the scalar power spectrum has an oscillatory behavior in the smaller value of k and a suppression of a 
large scale power. For the planar mode, fca <C k, the adiabaticity parameter vanishes in the earlier times the scalar 
power spectrum is well defined, but during the intermediate times it becomes greater than unity. Then, the effect 
of primordial anisotropy is enhanced and in the resultant spectrum the scale-dependence is unsuppressed. For the 
modes of k^ = 0, the adiabaticity parameter diverges and the WKB approximation is not well defined. But such a 
mode is not observable. 

In the case of the Bianchi IX model, by definition, the modes are discrete, but each mode exhibits a similar behavior 
to the corresponding mode in the case of Bianchi I model. For the moderate modes that satisfy K > 12io), where 
the quantum number K characterizes the angular momentum along the preferred direction k^/H in the Bianchi 
I model), in the resultant power spectrum the angular dependence is suppressed. Note that roughly speaking to 
characterizes the degree of the initial anisotropy, and i is related to the late time expansion rate of de Sitter {£ = 1/H). 
In contrast, for the planar modes K < £/{~12to), the effect of the primordial anisotropy is not suppressed. 

In summary, the anisotropy is not always suppressed in the present models. The anisotropy of the planar modes 
may leave non-negligible effects on the plane orthogonal to the preferred direction. They would be generic predictions 
from an initially anisotropic Universe in vacuum. There are important issues left for future studies. To obtain more 
reliable predictions, our method should be applied to solve the metric perturbations and in particular to find out the 
role of coupling of one of the tensor polarizations to the scalar mode. The evaluation of the bispectrum or trispectrum 
would also be important. 
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Appendix A: Evolution of Bianchi IX Universe 



Employing the Hamiltonian formalism, one can study the cosmology in the Bianchi universe in analogy with the 
classical mechanics. By parameterizing the scale factors of three independent spatial directions as 5"+'^++^'^-, 
ga+^+-\/3^_ and e"~^''+, the comoving evolution can be specified by two variables In analogy with the 

classical mechanism, the evolution of the Bianchi IX Universe can be described by the motion of a particle on the 
plane with the potential 

l/ix(/3+,/3-) = 26^^^+ cosh (4V3/3_) + e-^f^+ - 2e^'^+ - 46"^^+ cosh (2V3/3_). (Al) 

The potential is almost flat near the origin = 0,/3_ =0) but there is a steep exponential wall which forms an 
equilateral. Note that the corresponding potential in the Bianchi I model is vanishing, say Vi(/3+,/3_) = 0. The 
cosmic motion is almost confined inside this wall. There are three valleys into the potential along the straight lines, 
/3_ = > 0), /3_ = a/3/3+ and ^_ = -a/3/3+ (for /?+ < 0). Only along these valleys, the Bianchi IX Universe 
can evolve to have highly anisotropic geometry with large or These three valleys are equivalent under the 
exchange of label of the axes. 

The evolution of Taub-NUT de Sitter spacetime is initially along one of these valleys on the plane, 
starting from the infinity. In fact, one can read off the correspondence of the variables a + /?+ + \/3/3- = In a, 
a + /3-I- — = Ina, and thus a — 2/3+ = Inc, and one can easily find a = (1/3) In(a^c), /3+ = (1/3) In (a/c) and 

/3_ = 0. 
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